Observation of isotonic symmetry for enhanced quadrupole collectivity 
in neutron-rich 62,64,66pg isotopes at = 40 
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The transition rates for the 2^ states in 62,64,66pg ^gre studied using the Recoil Distance Doppler- 
Shift technique applied to projectile Coulomb excitation reactions. The deduced E2 strengths il- 
lustrate the enhanced collectivity of the neutron-rich Fe isotopes up to A'^ = 40. The results are 
interpreted by the generalized concept of valence proton symmetry which describes the evolution 
of nuclear structure around A'' = 40 as governed by the number of valence protons with respect to 
Z ~ 30. The deformation suggested by the experimental data is reproduced by state-of-the-art shell 
calculations with a new effective interaction developed for the fpgd valence space. 

PACS numbers: 21.10.Re, 21.10.Tg, 21.60.Cs, 27.50.+e 



Configuratioii mixing is one of the most important con- 
cepts of quantum mechanics. Orbital mixing of electron 
configurations in organic molecules provides extra stabil- 
ity to the system; quark mixing accounts for rare meson 
decays violating the CP symmetry. 

Configuration mixing in nuclear states provides a rich 
aspect of atomic nuclei, since it can considerably enhance 
amplitudes of collective excitation modes. Quadrupole 
deformation in nuclei is one such example and is known 
to evolve dramatically. Pronounced collectivity or defor- 
mation is normally observed in nuclei away from closed 
shells. It can be correlated with the dimension of the va- 
lence space, which is the product of valence proton and 
neutron numbers. This highlights a unique feature of 
atomic nuclei as a two fermionic quantum system. 

Recently, accumulated evidence has shown that magic 
numbers of nuclei can change far from the /3 stability 
line P, [2| • This has challenged our familiar picture of the 
shell closure, motivating the quest to understand the evo- 
lution of collectivity and mixing properties of associated 
configurations at extreme isospin. 

In the present work, we measured the mean lifetime (t) 
of the first 2+ (2^) states of the neutron-rich 62,64,66pg 
isotopes at and around the neutron number N = 40. As 
a direct measure of quadrupole collectivity, we deduced 
the reduced E2 transition probabilities i?(E2;2+^0"''), 
which are inversely proportional to r. Fe isotopes, span- 
ning from the most neutron-deficient two-proton emit- 
ter *^Fei9 0, via the most abundant isotope ^^Feao, to 



the most neutron-rich Fe48 hitherto observed [J|, pro- 
vide an intriguing testing ground for theoretical mod- 
els to comprehend the development and evolution of 
quadrupole collectivity over a wide range of neutron 
numbers. The i3(E2) values should follow an inverted 
parabola when plotted as a function of valence neutrons 
going from the closed shell at iV = 20 towards the next 
closed shell at = 28 and then again from A = 28 to 
A^ = 50. The magicity at A^ = 50 is still an open ques- 
tion in the regime of large neutron excess. A key issue is 
the degree of quadrupole collectivity manifested in ^^Fe4o 
at the A = 40 harmonic-oscillator shell closure. If the 
A^ = 40 shell closure is absent — as suggested by the low 
excitation energy £'(2+) of the 2^ state in ^^Fe |3| — one 
can expect the maximum collectivity at A^ sa 40, which 
corresponds to the middle of the magic numbers 28 and 
50. In fact, a recent study indicated a sudden increase of 
S(E2) from ^^Feae to ^'^Feag [H. However, this picture 
is in contrast to a doubly magic character of the neigh- 
boring A^ = 40 isotone ^*Ni4o, where a steep decrease of 
S(E2) has been observed from to ^^Ni 0|. 

In this Letter, we investigate the interplay between 
quadrupole collectivity in 62.64,66pg ^^j^^ ^j^g configuration 
mixing of valence nucleons in terms of the semi-empirical 
scheme of the valence proton symmetry Q as well as 
state-of-the-art shell model calculations. Based on new 
^^Fe data and additional ones that are more precise than 
those given in Q , we provide a new insight into the com- 
plex evolution of collectivity in the vicinity of A^ = 40. 
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The experiment was performed at National Supercon- 
ducting Cyclotron Laboratory (NSCL), Michigan State 
University. The lifetimes of the 2+ states in 62,64,66pg 
were measured by means of the Recoil Distance Doppler 
Shift (RDDS) technique applied to intermediate-energy 
Coulomb excitation. This technique has been demon- 
strated for medium- heavy nuclei with Z ?» 50 \ld^ . 
In this work, we extended the RDDS measurements to 
lighter Fe projectiles with Z = 26 by optimizing the 
target-degrader combination of the Koln/NSCL plunger 
device |9l4ll| . By making use of intense secondary beams 
provided by the Coupled Cyclotron Facility at NSCL, the 
Segmented Germanium Array (SeGA) for gamma-ray de- 
tection, and the S800 spectrograph for identification 
of reaction products, we realized the lifetime measure- 
ments of neutron-rich Fe isotopes up to = 40. 

Secondary beams of Fe were produced by fragmenta- 
tion of a primary ^^Ge beam at 130 AMeV incident on 
^Be targets. The A1900 fragment separator was used 
to purify the fragments. The resulting beam was sa85% 
^^Fe (at a typical rate of 3.6 x 10^ pps and an incident en- 
ergy of 97.8 ^MeV), «65% ^^Fe (6xl03 pps, 95.0 AMeV), 
and «25% ^^^Fe (1x10^ pps, 88.3 AMeV) for each setting. 

The beams were directed onto the plunger device, 
which was mounted at the target position of the S800 
spectrograph. The plunger target was a 300 /xm-thick 
Au foil, while a 300 fj,m (400 /xm)-thick Nb foil was used 
as a degrader for the 62,66pg (^64pg'j jneasurement. The 
beam velocities before and after the degrader were esti- 
mated to he v/c = 0.368, 0.322 for ^^Fe, v/c = 0.364, 
0.298 for 64Fe, and v/c = 0.346, 0.291 for ^^Fe, respec- 
tively, with typical velocity widths (in r.m.s.) of 
Mass and charge of the incoming beams were identified 
event- by- event using the time-of-flight measured as the 
timing difference between two plastic scintillators placed 
in the extended focal plane of the A1900 and in the object 
plane of the S800 analysis beam line. Scattered particles 
were identified by the detection system of the S800 (l2| . 

Doppler-shifted 7 rays were detected by SeGA coupled 
with the digital data acquisition system [l3 |. Two rings 
of 7 and 8 detectors were mounted at forward and back- 
ward angles of 30° and 140° relative to the beam axis. 
Data were taken at 5 — 7 different target-degrader dis- 
tances over the range from to 20 mm for different Fe 
isotopes. The data taken at the largest distances were 
used to examine the contributions from Coulomb excita- 
tion reactions in the degrader 0, [ll| . 

Lifetimes were determined by line-shape analysis in a 
similar procedure to those developed in Refs. 0,[lH. The 
best-fit results for each Fe isotope are shown in Fig. [T] 
The calculation takes into account the velocities of the 
incoming and outgoing beams, the energy losses through 
the target and degrader, and the detector geometry in- 
cluding effects due to the Lorentz boost. The energy and 
angular straggling of the reaction products and the en- 
ergy resolution of the 7-ray detectors are described as a 
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Figure 1. (color online) Typical results from the 7-ray line- 
shape fits to the data for the 2"''— transitions in 62,64,66pg^ 

single width parameter of a Gaussian distribution incor- 
porated in the lineshapes. We used four different width 
parameters for decays before and after the degrader in the 
spectra obtained with the forward and backward rings. 
The fits to the data were performed using variable param- 
eters which include the lifetime of the state of interest, 
the width parameters, a normalization factor, as well as 
the yield ratio between target and degrader. The large 
distance data were found to be consistent with degrader- 
to-target yield ratios of around 30-40 % for 62,64,66pg 

In the present analysis, special care was taken to select 
sensitive regions of the spectra which ensure the reliabil- 
ity of our analysis in determining the lifetime. To 
this end, after fixing the parameters that characterize 
the lineshape, we compared the simulated spectra to the 
measured ones by varying only the lifetimes. It became 
obvious that the most sensitive region corresponds to the 
7-ray energy range closest to the two peak centroids. It is 
worth noting that the peak amplitudes are most sensitive 
to the lifetime, whereas the value can also be affected 
by details of the lineshape. Therefore, we restricted the 
calculation to the optimum regions as indicated by 
the data with the error bars in Fig.[TJ We found a higher 
sensitivity to the lifetime for the spectra obtained from 
the forward-ring detectors. We thus first made a fit only 
to the forward-ring spectra to search for the minimum 
over a wide range of the lifetime. Later, we determined 
the lifetime from an overall fit around the obtained min- 
imum using the data taken with both rings. The results 
of T are summarized in Table HI together with the corre- 
sponding B(E2;2+— >0+) values. The final lifetime results 
are consistent with the first results obtained with the 
forward-ring data. For ^^^Fe, we used only forward-ring 
data due to a contaminant in the backward-ring spectra. 

The evolution of collectivity for even-even nuclei with 
Z = 24 - 36 and A^ = 22 - 58 is studied in Fig.[2]from the 



3 



Table I. Comparison between the experimental (Exp) and 
theoretical (SM) E{2+) and B{E2) values for 62,64,66pg rj,^^ 
lifetimes (texp) are the present results, while the E{2^) values 
are taken from Ref. [23|. Previous B(E2) data (Pre) @| are 
also compared. 
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Exp 


SM 


Exp Pro SM 


^^Fe 8.0(10) 877 


835 


198(25) 214(26) 270 


®*Fe 10.3(10) 746 


747 


344(33) 470ti"; 344 


^''Fe 39.4(40) 575 


570 


332(34) 421 



systematic behavior of the energy of the 2^ states (Fis 
(a)) and the B(E2;2+^0+) values (Fig. [U (b)) [ll-i? 
Only Ni isotopes reveal characteristics of shell closure at 
N = 28, 40 as inferred from large E{2+) and small B(E2) 
values. The = 40 shell closure vanishes for all the other 
even-even isotopes, where an enhanced collectivity is evi- 
dent. The present i3(E2) data confirm the onset of collec- 
tivity for ^^Fe and ^^Fe recently reported Q (see Table|I]) 
and represent the first evidence that the enhanced collec- 
tivity persists in ^^Fe at iV = 40. This tendency is con- 
sistent with the other neighboring isotopes with Z ^ 28. 
A large deformation parameter of /32 ~ 0.28 is calculated 
from the measured B(E2) of '^'^'^^Fe. 

Interestingly, one can see in Fig. [2] (a) no- 
table similarities in the £'(2+) behavior for 
the pairs of Fe {Z = 26)/Se (Z = 34), and 
Zn {Z = 30) /Ge {Z = 32)^. An incipient similarity in 
E{2+) is also suggested for the Cr (Z = 24)/Kr {Z = 36) 
pair. The enhancement of collectivity is also noticeable 
from the decrease of i?(2+) as one moves away from the 
proton number Z w 30. It starts at = 36 for Cr/Kr, 
iV = 38 for Fe/Se aird = 40 for Zn/Ge, respectively. 
These striking similarities suggest an isotonic symmetry 
in the developmeirt of collectivity for each of these pairs. 

To further investigate the underlying symmetry in- 
ferred from the £^(2+) behavior, we turn to the valence 
proton symmetry (VPS) [1, i]. The VPS scheme cor- 
relates a pair of isotones that have the same number of 
valence proton holes and particles with respect to closed 
shells. The same valence space for correlated nuclei re- 
sults in a similar degree of collectivity and its evolution, 
as shown, for example, in the region of 50 < < 82 and 
.Z « 50 [9| . In this work, we generalized the VPS concept 
since the symmetry appears with respect to the proton 
number Z = 30 instead of the canonical magic number 
Z = 28. Here, we consider proton holes in the I/7/2 or- 
bital below Z = 28 and protons in the I/5/2 orbital above 
the Z = 32 subshell closure. This is justified by the fact 
that the 2^3/2 orbital is isolated between the I/7/2 and 
1/5/2 orbitals. The 2^3/2 orbital, if it is fully occupied, 
has no impact on the evolution of collectivity, and hence, 
for the VPS scheme, the effect of a fully empty or occu- 
pied 2p3/2 orbital is identical. Consequently, Se (Kr) can 




Figure 2. (color online) (a) Energies -5(2"*") and 

(b) _B(E2;2+^0+) values for even-even nuclei with Z = 24 - 
36 and N = 24 - 56. The stars indicate the present data 
of 62.64.66pg Tj^g 5(E2) data of Kr in (b) are divided by 2. 
In (c), the _B(E2) data of the Fe and Cr isotopes are com- 
pared to the scaled _B(E2) values for the Se and Kr isotopes, 
respectively, as well as the shell model predictions (SM) for 



62,64,66,68 



Fe and Cr isotopes. 



be a valence partirer of Fe (Cr) . 

To test the modified scheme of the VPS, in Fig. [5] (c) 
we compare the existing i?(E2) data for the Fe/Se and 
Cr/Kr pairs, including the present data of 62,64,66pg^ rpj^g 
charge and mass differences are considered iir the compar- 
ison by scaling i?(E2) data for the heavier isotopes (Se 
and Kr) with the factor of S" = {ZL/ZH)^^iAL/AH), 
where Zh and Ah {Zl and A^) denote the atomic and 
mass numbers of heavier (lighter) pair partners, respec- 
tively 0. While the data are still scarce for Cr/Kr, the 
present data of 62,64,66pg clearly agree in magnitude with 
the overall trend of the Se isotopes, proposing the validity 
of the generalized VPS scheme. Consequently, one can 
expect large collectivity for further neutron-rich Fe and 
Cr isotopes beyond A^ = 40, as illustrated by the scaled 
B{E2) data for the Se and Kr isotopes. 

The symmetric evolution of collectivity above and be- 
low the proton number Z = 30 can be understood in 
a shell- model picture [23|, in terms of an energy de- 
crease of the neutron shell gap at A^ = 40 caused by 
the proton-neutron monopole tensor interaction. Nuclei 
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above Z = 30 become more deformed when approaching 
N = Z = 40 ^°Zr. Filhng protons in the I/5/2 orbital 
lowers the neutron lgg/2 and 2^5/2 orbitals. Vice versa, 
for the light isotones with Z < 26 that have proton holes 
in the I/7/2 orbital, the attractive interaction with the 
neutron I/5/2 orbital and the repulsive interaction with 
the 1(79/2 and 2^5/2 orbitals get weaker. In both cases, 
the shell gap at iV = 40 reduces consistently, enhancing 
quadrupole collectivity. A similar effect is produced on 
the Z = 28 gap. This is consistent with the ground-state 
spin parities of the = 39 isotones [l^ ( ^^Sr : 5/2"*", 
75Kr : 5/2-1", ^SSe : 9/2+, ^^Ge : 1/2", ^^Zn : 1/2", 
^^Ni : (l/2)~ ), where the intruder 1^9/2 (or 2^5/2) con- 
tribution is evident only for the heavier isotones away 
from Z = 30 that have positive-parity ground states. 

To allow a quantitative discussion of the intruder neu- 
tron configurations for the lighter valence partners, we 
have performed shell model calculations [28| and com- 
pare them with the present B(E2) data for 62,64,66pg^ 
The natural choice for the model space is the fp shell 
(20 < Z < 40) for protons and the 2p3/2l/5/22pi/2lff9/2 
space for neutrons (28 < N < 50). However, shell model 
calculations in this model space fail to fit the ^®Fe 2^ 
state [l^l, while they can reproduce the level schemes of 
62-64pg j-ather well. To reproduce the quadrupole col- 
lectivity in this mass region, the inclusion of the neutron 
2(^5/2 orbital is necessary 0,[2l|- Recently, a new effective 
interaction for this large model space (LNPS) has been 
developed that can reproduce the different phenomena 
suggested by the available data in this mass region (ist . 
The LNPS interaction has been obtained by adapting 
a realistic CD-Bonn potential to this model space after 
many-body perturbation techniques j29| and monopole 
modifications. The results for the excitation energy and 
B(E2) values obtained with standard effective charges 
(e„ — 0.5e,ep = 1.5e) for the Fe isotopes of interest are 
reported in Table HI 

The excitation energies are very well reproduced by 
the calculations in all cases. The calculated i?(E2) of 
^^Fe coincides with the data, while overestimated values 
are obtained for ^^'^^Fe. The theoretical i3(E2) values 
increase linearly with the neutron number but the exper- 
imental values seem to suggest a slightly fiatter behavior. 

The analysis of the ground state wave functions shows 
that in ''^Fe, the probability of excitations to the gd or- 
bitals is about 60%, increasing to - 90% in S'^Fe and to 
almost 99% in ^^Fe. This is consistent with the devel- 
opment of a new "island of inversion" at iV = 40 in this 
mass region. The calculated B(E2) values of 62.64,66,68pg 
and ^°Cr are plotted in Fig.[2](c). The trend is similar to 
that expected for these isotopes from the VPS scheme, 
suggesting that a region of enhanced collectivity extends 
to neutron-rich Fe and Cr isotopes beyond N = 40. 

In summary, we measured the lifetimes of the 2f states 
of 62,64,66pg^ quantifying for the first time an enhanced 
collectivity of the Fe isotopes at TV = 40. The observed 



trends in £^(2+) and i?(E2) are interpreted by the va- 
lence proton symmetry with respect to the proton num- 
ber Z « 30, suggesting similarities in the role played by 
active protons or holes in the I/5/2 and I/7/2 orbitals, 
respectively. The shell model calculations with the new 
effective LNPS interaction describe the present B(E2) 
data well. The calculations point to the important role 
of the neutron intruder configurations where excitations 
to the Igg ji and 2d^ /2 orbitals across the iV = 40 subshell 
gap are favored due to the shell migration far from stabil- 
ity. These results suggest that, owing to the collaborative 
interplay between the valence proton and intruder neu- 
tron configurations, a new deformation region appears in 
the neutron-rich Fe isotopes at and beyond iV = 40. 
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